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ABSTRACT 

Oh; 

Q ■ We have obtained deep surface photometry of the edge-on spiral galaxy 

+i . NGC 4244. Our data reliably reach 27.5 R magnitude arcsec"^, a significant 

d ■ improvement on our earlier deep CCD surface photometry of other galaxies. 

NGC 4244 is a nearby Scd galaxy whose total luminosity is approximately one 
S^ . magnitude fainter than the peak of the Sc luminosity function. We find that it 

has a simple structure: a single exponential disk, with a scale height hz = 246 
± 2 pc, a scale length hji = 1.84 ± 0.02 kpc and a disk cutoff at a radius Rmax 
= 10.0 kpc (5.4 scale lengths). We confirm a strong cutoff in the stellar disk at 
Rmax, which happens over only 1 kpc. We do not see any statistically significant 
evidence for disk flaring with radius. Unlike the more luminous Sc galaxies 
NGC 5907 and M 33, NGC 4244 does not show any evidence for a second 
component, such as a thick disk or halo, at yU/j < 27.5 magnitude arcsec"^. 

Keywords: galaxies: individual (NGC 4244) — galaxies: photometry — 
galaxies: spiral 
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1. INTRODUCTION 



In order to better characterize the distribution of luminosity in the stellar components 
of disk galaxies, we are conducting a detailed study of edge-on disk galaxies with a range of 
bulge-to-disk ratio using deep CCD surface photometry. Our galaxies range in Hubble type 
from SO to Scd, and in bulge-to-disk-ratio from 20 (NGC 3115) to 0.03 (NGC 4244, this 
work & [Baggett et al. 1998| ). We will re-examine the structure of their bright disks and 



bulges, and also study, for the first time, the structure of faint outer components such as 
thick disks and stellar halos. 

These faint outer components are important because they may trace the history of the 
galaxy, either because they contain the oldest stars ( |li)ggen et al. 196^ , [h^dvardsson et al. 



|1993|) or because they were formed via accretions of smaller galaxies ( pearle fc Zinn 1978|, 



Davis et al. 19851 IQuinn fc Goodman 1986|) . In addition, if they are made of old stars, they 



should provide connections to high-z galaxy observations, both direct (|Steidel et al. 199^ ) 
and indirect via absorption line studies ( [Wolfe 1990|) . 

The simplest galaxies structurally are those of late type - Sc and later. To first order, 
they comprise simply a disk with exponential surface brightness falloff ( |I''reeman 1970] , van 



der Kruit & Searle 1981). We have already studied a large Sc galaxy, NGC 5907 (Morrison 
et al. 1994). Our analysis confirmed previous findings on its thin disk (van der Kruit & 
Searle 1981), but we found an unusual second component, whose luminosity falls off very 
slowly with radius (see [Sackett et al. 199^ , [Monkiewicz et al. 1999|) . Since the galaxy 



luminosity function (LF) for Sc galaxies is very broad, and peaked at total blue luminosity 
iMB,tot) of -19.5 (IBinggeli et al. 1988|) , making NGC 5907 iMB,tot = -20.55, [fully 1988|) 



an unusually luminous Sc galaxy, we have also studied NGC 4244, a smaller late-type spiral 
whose total magnitude of -18.4 B mag ( piling 1996a| ) is almost a magnitude below the 



broad peak of the field Sc LF. NGC 4244 has Hubble type Scd, is located 3.6 Mpc away 
(Tully-Fisher, [Aaronson et al. 198tj| ), and has a fiat rotation curve with Vrot = 100 km/s 
( piling 1996a|) . Beyond the optical disk, NGC 4244 has a declining rotation curve ( piling| 



[1996 



a 



In this paper we present deep surface photometry of NGC 4244. Our data reach 
reliably to fiR = 27.5 magnitude arcsec"^, significantly fainter than our previous published 
data on NGC 5907. Questions that we will address using this extremely deep data are: 

• do small late-type disk galaxies such as NGC 4244 have similar properties to large 
ones such as NGC 5907? 

• what are the properties of NGC 4244's luminous thin disk? Does the stellar disk 
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have a distinct edge? What is its vertical scale height? Does this scale height remain 
constant with radius? 

• does NGC 4244 show fainter components like a thick disk or halo, or like NGC 5907's 
unusual second component? 

NGC 4244 is less massive than a typical Sc: the total mass of NGC 4244 is 9.77 x 
IO^Mq ( [Tully 19"88|) , whereas the average total mass of an Sc is ~ IO^^Mq ( [Roberts"^ 



Haynes 1994 ). However the H I mass fraction Mni/Mtot ( piling 1996a ) is fairly typical 



for an Sc galaxy: 0.14. The mean Mm/Mtot for a typical Sc (from [Roberts fc Haynes] 



T^MD is 0.1. The mass fraction of Ha {MujUtot) is low in NGC 4244, relative to other 



galaxies of the same mass. In his CO study, [Sage (1993)| found < Mh2/Mt >= 0.0013 



which placed NGC 4244 in his mass class HI, for which < Mh2 > = 8.4 x 10'' Mq or 
< Mh2/Mt >~ 0.008. The current star formation rate in NGC 4244 is also low, as can be 
seen from its low X-ray, radio continuum, and IRAS emission (cf. [Oiling 1996aD . 



The Hubble type of a galaxy is known to correlate roughly with its environment in the 
sense that late- type galaxies are found in more isolated areas ( pressler 1980| ). In addition, 
we might expect to see a correlation between the existence of thick disks and galaxy 
environment if the accretion theory of thick disk formation is correct. NGC 4244 lies in 
the Coma-Centaurus spur, and is thought to be a member of the CVn I ( [de Vaucouleure 



1975 ) or B4 ( Kraan-Kortweg fc Tammann 1979 ) group. Its small distance makes group 



membership very difficult to decide, since group peculiar velocities are of the same size as 

the Hubble flow. We have plotted in Figure [| supergalactic coordinates of galaxies in the 

B4 and nearby B5 groups of Kraan-Kortweg & Tammann - it is not clear whether these two 

groups are separate entities. We can only conclude from these data that NGC 4244 lives in 

a relatively isolated environment, with only three galaxies of similar or greater brightness 

within ~ 2 Mpc (NGC 4449, NGC 4736, NGC 4826), the brightest of which (NGC 4736) 

has Ms^tot = -19.37 mag, 1.7 mag brighter than NGC 4244. 

I ' . „ ^1 

van der Kruit & Searle (1981a) used photographic surface photometry of NGC 4244 

to study its structure, flnding the disk well modeled by a radial exponential function and 

a vertical sech^ function, with a sharp cutoff at 5.3 radial scale lengths. They found no 

evidence for a thick disk, in contrast to their studies of disk galaxies with larger B/D 

ratios where thick disks were seen. However, thick disks are very faint, compared to both 

the night sky and the light of the younger, thin disks, and, thus, detection and accurate 

determinations of the properties of thick disks from photographic work are difficult. The 

improved precision of CCDs is needed, and our faint CCD photometry of NGC 4244 will 

put much stronger limits on the existence of a thick disk in this galaxy. 
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The flaring of a galaxy's stellar disk gives constraints on the balance between the 
gravitational potential the disk experiences and heating mechanisms the stars experience, 
van der Kruit & Searle (1981a,b,1982a,b) found that stellar disks in the galaxies they 
studied had constant scale heights with radius, although this has recently been challenged 
(de Grijs & Peltier 1997). NGC 4244 is a good galaxy in which to study disk flaring, as it 
has no large companions and no apparent thick disk (cf. van der Kruit & Searle 1981a). 
We will discuss the strong limits our data place on any variation of disk scale height with 
radius. 



2. OBSERVATIONS & REDUCTIONS 

The NGC 4244 data were obtained over a five night run in March 1997 and a six night 
run in April 1997 using the 24"/36" Burrell Schmidt telescope^ at Kitt Peak National 
Observatory]^. The images were taken using a 2048x2048 frontside-illuminated Tektronix 
CCD (S2KA) mounted at Newtonian focus. With this setup, the field of view is 69 square 
arcmin, with each 21 /im pixel imaging 2. 03 of sky. The gain was set to 2.5 e~/ADU and 
the readout noise was 3 e~ (1.2 ADU). All exposures were made through a Harris R filter. 



Observations were performed in the manner described by [Morrison et al. (1997) . Half 



of the time was used observing the galaxy and the other half was used to obtain dark 
sky flats. The dark sky flat images were taken at approximately the same hour angle and 
declination as the galaxy images. Over the course of the two observing runs, 66 galaxy 
images and 70 dark sky flats were obtained. For all of the galaxy and sky images, the 
exposure time was 600 s. Conditions were photometric throughout the observing runs. 

Overscan removal and bias subtraction were done in the standard manner. Then 
"master" sky flats were constructed from the dark sky images from each run. The individual 
sky flats were dithered by ~ 5°, so that wings from bright stars in the field would not 
appear in the final, combined image. To construct the master sky flat, we used an iterative 
procedure. We scaled the individual sky images by their mode. Prescaling was important 
because, even during photometric conditions, the sky brightness varied by ~0.3 mag. This 
was likely due to the long tail to high wavelength of the Harris (Cousins) R filter, which 



■^Observations made with the Burrell Schmidt of the Warner and Swasey Observatory, Case Western 
Reserve University. 

''Kitt Peak National Observatory is a division of the National Optical Astronomy Observatories (NOAO), 
which are operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under 
cooperative agreement with the National Science Foundation 
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includes some of the strong airglow lines, combined with the change in sky brightness with 

zenith distance. Then, we combined the individual, scaled, sky images, using IRAF's ccdclip 

algorithm, which uses CCD properties to estimate expected variations in data values and 

removed pixels which differed from the median by more than 2cr. Each of the individual 

sky frames was divided by the flatfield frame to reduce the width of the distribution of sky 

values and so make rejection of outliers due to faint stars and stellar wings more accurate. 

A plane was fit to each of the individual flat-fielded sky frames and subtracted from the sky 

frame. The modes were recalculated, and the entire procedure repeated. This procedure 

was repeated until no improvement was seen in the master sky flat (ten times in this case). 

// 
The 66 galaxy images were dithered by ~ 400 , much less than for the sky fiats. The 

individual galaxy images were fiatfielded using the master sky fiat and registered. We 

combined the individual galaxy images, using a median. In Figure |^ we show our final 

galaxy image, after subtracting a best fit plane sky and correcting to airmass 1. 

One Landolt standard star field ([Landolt 1992| ) was observed each night, giving us 26 
well-exposed standard stars with a range in color and airmass. A photometric zero point of 
R = 22.23 ± 0.03 magnitude arcsec"^ corresponding to 1 ADU/sec/pixel was determined. 
For a 600 s exposure, this yields R = 29.18 ± 0.03 magnitude arcsec"^ corresponding to 1 
ADU/exposure/pixel. As our exposures were only taken in one filter, we cannot add a color 
term to our calibration but estimate its magnitude as less than 0.10 mag over the color 
range of our standards (see Figure ^ B — V = —0.06 to 1.66 mag, V — R = 0.05 to 0.8). In 
fact, over the expected color range of the relatively old populations we will be studying, the 
color term is well below the uncertainty in the zero point. 

The average sky brightness, 1998.0 ADU / pixel (see Section 3.2), corresponds to Rgky 
= 20.9 magnitude arcsec"^. 



3. ANALYSIS & RESULTS 

For simplification in subsequent analyses, the galaxy has been rotated counter-clockwise 
by 42°. This agrees well with Olling's (1996a) determination of —48° for the position angle 
of NGC 4244. 
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3.1. Masking 



In order to reach the faintest possible surface brightness level in the galaxy, we masked 
out the foreground stars and background galaxies in our field. To characterize the stars we 
used DAOphot ( Stetson 1987|) to determine the positions and magnitudes of all of the stars 



within the frame. Using 48 bright, but unsaturated, stars, we determined the point spread 
function (PSF) for stars in the frame. We followed the method described in Stetson (1987), 
iterating until a clean subtraction within a 15x15 pixel box centered on the PSF star was 
obtained. Our aim was not to determine accurate photometry of all of the stars within the 
frame, but to mask the stars out to a constant level in surface brightness. This means that 
we paid particular attention to the outer wings of our PSF. We used this PSF to construct 
a radial profile and masked all pixels for which a foreground star contributed at least 1 
ADU. The stellar radial profile was well determined out to the PSF radius, 15 pixels (30"). 
The PSF is axisymmetric and varies along the chip. The variation is most significant for the 
central regions of the stars and does not strongly affect the outer wings of the stars, which 
are formed by processes such as scattering in the atmosphere and telescope. We describe 
the PSF beyond 15 pixels in Section 3.4. 

Background galaxies and undetected stars were masked by hand using the IRAF task 
imedit. Special care was taken with stars near the galaxy, where the galaxy's luminosity 
gradient interferes with DAOphot's sky fitting algorithm. 

NGC 4244 has a small dust lane and several visible H II regions. To mask these, a dust 
mask 20 pixels wide, running along the major axis, was created. In Figure ^ we show the 
galaxy with our final mask superimposed. 

Using this mask, the data were binned up into regions whose size varied from a 
resolution element near the galaxy's center to 50 x 50 pixels in the regions far from its 
center. This procedure used a robust averaging scheme to compensate for any unmasked 
stars, galaxies, etc. and is described in detail in [Morrison et al. (1994) . 



3.2. Sky Subtraction 

Accurate sky subtraction is crucial for the faint isophotes. To determine the sky level, 
we used the masked image and made a histogram of sky values in bins well away from the 
galaxy, using bins with more than 100 pixels after masking. There are 310 such bins. Using 
a histogram of the sky values within these bins, we determined the sky value to be 1998.0 
±0.8 ADU. In Figure ^ we show this histogram. The width of the histogram provides a 
measure of the contribution of unmasked stars to the sky level. 
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3.3. Limits to Our Precision 

Our photometry reaches rehably to 27.5 magnitude arcsec"^. The error model is 
described in Appendix A. To illustrate, we work through the errors in a 13 x 3 pixel bin, 
centered 12 pixels (425 pc) from the major axis and 150 pixels (5.3 kpc) from the minor 
axis. In this bin, the mean number of counts is C = 2498 ADU per pixel, with the mean 
number of counts from the galaxy alone y = 500 ADU per pixel. The errors are summarized 
in Table |I|. The dominant source of error in the thin disk region of the galaxy is the error 
arising from surface brightness fluctuations within NGC 4244, because NGC 4244 is so 
close to the Milky Way (its distance is 3.6 Mpc; note that it is marginally resolved even 
from the ground.) Surface brightness fluctuations dominate the formal errors wherever 
'g/upixeis > 0.40 (i.e. up to 1.6 kpc above the plane). In the faintest region of the galaxy 
(above 1.6 kpc from the plane), the main sources of error are large-scale fiat-fielding errors 
and sky variations. 



3.4. Scattered Light? 

To see how much scattered light from the central part of the galaxy contributes at 
high-z, we constructed a stellar profile out to 150 pixels (5.3 kpc). In the center, we 
used the stellar profile constructed using DAOphot (Section |3]5). For the outer region, 
we used two saturated stars. Using stars whose centers were saturated means that the 
signal-to-noise in the outer wings is higher than in unsaturated stars. We masked the stars 
in the wings of the saturated stars, as well as diffraction spikes around the saturated stars. 
The unmasked pixels were averaged in radial annuli. The profiles from the two saturated 
stars were averaged and joined to the DAOphot PSF profile. The combined profile was 
normalized to the surface brightness of the central bright region of NGC 4244. The final 
profile is shown in Figure ^. At a distance of 2 kpc (56 pixels) from galaxy center, the 
scattered light from the bulge is ~ 28.5 magnitudes arcsec"^, below our quoted photometric 
limit of hr = 27.5 magnitudes arcsec"^. Thus, it is not a significant effect. 



3.5. 2-D Fits 

We fit a single exponential disk to the data: 

Pd^.fc = i^oe-^/"«e-l^l/'^^ {R<Rma.y, (1) 

= iR>R„,a..) (2) 



(here Lq is the central density of the disk in Lq pc ^). The two-dimensional modeling is 
described fully in [Morrison et al. (1994)| . We start with a full three-dimensional model, 



and then for each iteration of the fit, we numerically integrate along the line of sight and 
produce a two-dimensional model, which we compare with the observed surface brightness. 

We start out with eight parameters: the scale length and scale height hji and hz, 
inclination i, the coordinates of the galaxy center xc and yc, the radial extent of the galaxy 
Rmax, the central surface brightness Ic, and the sky value. Because of high correlations 
among the parameters, eight free parameters are extremely difficult to fit accurately, but, 
fortunately, some of the parameters can be fixed using other information. The sky value is 
determined by averaging the values in any bin with more than 100 pixels, away from the 
galaxy, as described in Section 3.2. The xc and yc position of the center of the galaxy and 
Rmax are determined from a contour plot of the galaxy and by making one-dimensional 
profiles across the minor axis and parallel to the major axis. The inclination i = 84.5° was 
adopted from OUing's (1996a) kinematic measurement. Only three parameters then need to 
be determined from the 2-D fitting of a single exponential disk: hji, hz, and central surface 
brightness. 

How well do 2-D fits work? Because we are using a Levenberg-Marquardt algorithm 
( press et al. 1992|) , which is an iterative nonlinear leastsquares algorithm, we need to 
provide starting values. We investigated the sensitivity of the algorithm to the initial 
estimates used by mapping out the value of x^ over a large grid of starting values and 
noting where the algorithm converged. We find that if all three parameters, hn, hz, and 
Ic, are estimated, then the initial estimates need to be within 20% of the correct values 
for the model fit to converge. If one of the three parameters can be estimated from other 
information, such as hz from the 1-D fits, then initial estimates of the other two can be as 
much as 50% wrong and the model will still converge successfully. A starting grid of x^ 
can be used as a "Levenberg-Marquardt-by-eye" to locate the "trench" , that is, the starting 
parameter values for which the model fit will converge. Once the trench has been found, all 
three parameters are allowed to vary. 

As illustration, in Figure |^, we show a contour plot of x^ for fixed hz = Q.5 pixels (233 
pc). The contours are logarithmic; steep slopes are even steeper than they appear. For all 
starting values of hR and Ic shown on the plot, except upper left {h^ ^ 1000, Jc ~ HO), 
lower left {hji i=a 1000, Ic ~ 30), and lower right {hji ^ 3000, Ic ~ 30) corners, the starting 
model will converge. However, for hz values which are more than ^ 20% from the true hz, 
the initial values of hji and Ic need to be within !^50% of the true values. 

Another concern is the reliability of the estimates once the algorithm has converged. 
Morrison et al. (1994)| showed that estimates based on the covariance matrix of the fit 
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can underestimate the true values severely. Thus, to determine error estimates on the 
derived disk parameters hji, hz, and Ic, we used bootstrapping with resampling ( |Et'ron 



rgg^i Chapter 15 of [Press et al. T^ , [tili'ron fc Tibshirani 1993| ). We ran our 2-D model 



fitting program 1000 times. Each time, the data to be fit were chosen by random sampling 
with replacement from the actual data. The sets of resulting disk parameters were roughly 
Gaussian distributions, centered on the result of the run using all of the data. The 
errors estimated from the widths of bootstrapping analysis are more reasonable estimates 
of our true errors, as the bootstrapping errors include not only the formal errors from 
known contributions but also a more realistic estimate of goodness of fit and systematics 
(unmasked foreground stars, dust, and small, faint, background galaxies). In Figure |, we 
show a histogram of hz values estimated by our bootstrapping code. 

In Table ^ we show the results of our 2-D fitting with the bootstrapping errors, as well 
as the correlation coefficients between all the variables. To determine all of the correlation 
coefficients, all eight variables were started at their correct values and allowed to vary. The 
important correlation coefficients are those between the first three variables, the scale length 
hji, the scale height hz, and the central surface brightness /^oQ But the strong correlations 
between some of the other variables, such as between hz and inclination, or Rmax and h^, 
demonstrate why it is important to use other methods to estimate some of the parameters. 

We estimated Rmax directly from the image to reduce the effects of these correlations. 
Unfortunately, the galaxy is not axisymmetric: Rmax on the southwest side of NGC 4244 
(left side of Figure 4) is 8% larger than on the northeast side: 10.5 kpc vs 9.7 kpc. We chose 
an intermediate value of 10.0 kpc. Using the northeast side (right side of Figure 4) only in 
our estimation, with Rmax = 9.7 kpc, results in an estimate of hji which is 3% larger. The 
estimates of hz and /iq change an insignificant amount. The bright stars and H II regions 
at the left edge of NGC 4244 make it difficult to model the left side only. 



3.6. Vertical Profiles 

In Figure ^ we show vertical profiles parallel to the minor axis. Our best fit model 
of a single exponential disk with scale height of 246 pc is shown with the data. A single 
disk fit is a very good match to the data down to 27.5 magnitude arcsec"^ and there is 



^Note that the central surface brightness quoted here is for the observed inchnation of 84. °5. For 
comparison to other galaxies, the face-on and edge-on central surface brightnesses for this model are also 
quoted. 
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little evidence for disk flaring, i.e. changing scale height as a function of radius. The only 
profile that appears to flare is the one at ± 8.5 kpc. However, the warp in NGC 4244 starts 
at ^ 8.0 kpc and we attribute the apparent flaring to the warp. In the galaxy center, 
we may have a small contribution to the light from the bulge, which may lead to a slight 
underestimate of hz there, due to the steeper bulge profile. 

To test the constancy of hz with radius, we split the galaxy into inner and outer 
regions, R < 4 kpc and 4 < i? < 8 kpc (the region outside R =8 kpc was not included to 
minimize the effect of the small warp), and analyzed the two regions separately. For the 
inner region, we find hz = 246 ± 2 pc and, for the outer region, we find hz = 259 ± 8 pc. 
The difference between the inner and outer regions is not significant at the 2a level. This 
lack of evidence for flaring will be discussed further in Section 4.3. 



In Figure |10| and [Tl|, we show profiles parallel to the major axis, at ± 760 pc and at ± 
955 pc from the major axis. The horizontal profile at z = ± 760 pc is the closest distance 
to the major where a complete profile can be drawn. This is because of the large number of 
pixels masked due to dust, H II regions, and foreground stars at lower z. In Figure 0, the 
region further out than ~ 7 kpc from the minor axis shows more light than predicted by a 
single exponential disk fit, as if either the outer regions have a larger scale length h^ or the 
inclination changes slightly with radius, becoming closer to edge-on for large R. OUing's 
(1996a) H I data also suggest that the inclination increases slightly with radius. 

The profile at z = ± 955 pc shows the influence of the warp more clearly. Oiling 
(1996a) found evidence for a warped H I disk. He found that the warp starts close to the 
edge of the optical disk {R = 8.5 kpc), as is typical of galaxy warps ( |Briggs 1990| ). The H I 



warp extends below the major axis on the left side of the disk and above on the right side in 
our Figure 3. piling (1998)| finds that its line of nodes is close to the line of sight. In Figure 



[TT], the warp is apparent for the entire outer half of the disk, starting at 5 kpc from the 
minor axis, where data from the four quadrants of our image are quite different. In other 
words, in our faint optical data we see more details of the warped outer disk superimposed 
on the line of sight than is seen in the H I data, despite its high quality. (We do not show 
data past 7 kpc on the left side because of the presence of a large region of current star 
formation and three bright stars.) 

We also see evidence for a small warp in the outer parts of the stellar disk in our 
vertical profiles. The line of centers in the vertical profiles changes by ^2.5 pixels (89 pc) 
between the center and the left edge of the disk. We corrected for this effect before making 
Figure 9. This is equivalent to a change of 0.°5. The optical warp we find is in the same 
direction as the H I warp. The optical warp in NGC 4244 is much smaller than in NGC 5907 
(Morrison et al. 1994) or NGC 4565, which permits accurate study of its outer stellar disk. 
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4. DISCUSSION 



4.1. Comparison to Previous Surface Photometry 

I . ^ n 

van der Kruit & Searle (1981a) used photographic photometry to obtain thin disk 
parameters for NGC 4244 and, recently, Olhng (1996a) reanalyzed their photometry to 
rederive a disk scale height and disk scale length. We confirm the existence of a disk 
cutoff, first noted by [van der Kruit fc Searle (1981a)| , with our deeper data. [Kodaira "^ 



Yamashita (1996")| used J and K' photometry to estimate thin disk parameters. We list 
our disk parameters with those of the other groups in Table ^ including the orientation 
("as-observed", face-on, and edge-on), [van der Kruit fc Searle (1981a)| adopted a distance 
of 5 Mpc, rather than 3.6 Mpc used by the other studies, and, in Table ^, we have scaled 
their scale height and scale length to the distance adopted by the other studies. To compare 
their Bj data to our R photometry, we assume B — Bj = 0.2 and B — R = 1.5, colors 
typical for an old stellar population ( Bottema 1989| ). van der Kruit fc Searle (1981a) fit 



their data to a sech^ model, rather than an exponential light distribution and we include an 
estimate of the central surface brightness from an exponential model: /i (exponential) = fi 
(sech^) - 0.75 mag. 

IF we compare our R photometry to van der Kruit & Searle's Bj photometry, we find 
a blue color {B — R = 0.2 mag, rather than 1.5 mag) for NGC 4244. However, it is very 
difficult to determine accurate zero points from photographic photometry; B CCD data are 
needed to determine the galaxy's color. 



Our disk scale height is 15% larger than found by [van der Kruit fc Searle (1981a)| and 



by Oiling (1996a), but our scale length is the same. Our disk parameters do not agree well 
with those found by [Kodaira fc Yamashita (1996)| , however, their data are limited to the 



inner 180" (3.2 kpc), making the scale length more difficult to determine. In summary, our 
results agree well with previous optical studies. 

In Table H, we compare our disk parameters for NGC 4244 to those of NGC 5907 
(Morrison et al. 1994) and the Milky Way (cf Morrison et al. 1997). NGC 4244 has a 
smaller disk scale height than the other two galaxies and a fainter central surface brightness. 
The ratio of disk cutoff to radial scale length is similar to that of the Milky Way. 



4.2. Disk radial cutoff 



van der Kruit & Searle (1981a) showed that stellar disks do not stretch out indefinitely, 
showing a sharp cutoff at ~5 radial scale lengths, well inside the H I disk. Our deep 
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surface photometry of this almost unwarped stellar disk gives a strong test for disk cutoffs. 
Figure [1^ shows profiles just inside Rmax (< 10 kpc) and outside it (11.3 kpc). All the 
points outside Rmax have ^r fainter than 27 magnitude arcsec"^. The cutoff in stellar light 
happens relatively rapidly, over a distance of 1 kpc, and no faint outer component is visible 
down to iiR= 27.5 magnitude arcsec"^. The existence of this stellar cutoff is consistent with 
the detection by [Ferguson et al. (1998)| of faint H II regions "at distance twice the optical 
radius" of the nearby face-on galaxies they studied. The "optical radius" (corresponding to 
a face-on surface brightness ^b = 25 magnitude arcsec"^) is an outmoded concept based on 
what could be detected on photographic (POSSI) sky survey plates; it is now possible to 
see much fainter with CCD detectors. (For edge-on galaxies the situation is, of course, even 
easier.) Assuming hfj/hz = 10 to convert face-on to edge-on surface brightness, fiB = 25 
corresponds to /i/j ~ 21. Figure ^ shows that much of the stellar luminosity of NGC 4244 
comes from outside the "optical radius" defined in this way. 



4.3. Constant Scale Height 

What determines the thickness of a stellar disk? Since its brightness varies exponentially 
with radius and the contribution of dark matter should not be strong in the inner galaxy, 
the weaker gravitational pull from the outer disk means that we would expect a component 
whose vertical energy does not vary with radius to flare towards the edges. This can be 
seen clearly for NGC 4244's H I in Fig. 13 of Oiling (1996a), the scale height of H I 
increases from f» 400 pc at a radius of 5 kpc to 1.5 kpc at 13 kpc. However, the factors 
which determine the scale height of a stellar disk are more complex: although stars are 
born out of ISM with low velocity dispersion which is constant with radius (Pronfman et 



al 1988| , [Wouterloot et al. 1990| , [Scoville et al. 1993| , [Combes fc Becquaert 1997]) , various 



heating mechanisms change the velocity dispersion and, thus, the scale height. The velocity 
dispersion of stars in the solar neighborhood increases with age ( pahreiss fc Wielen 1983[ ). 



This is thought to be due to heating by spiral arms and giant molecular clouds, with GMCs 
being particularly important for vertical heating ([Jenkins fc Binney 1990[) . van der Kruit & 
Searle found that the scale height of the stellar disks in the galaxies they studied (1981a,b, 
1982a,b) remained constant with radius, despite the variation in CO surface density and, 
so presumably, numbers of GMCs with radius. 



Recently, [de Grijs fc Peltier (1997)] claim evidence for disk flaring in their study of 48 



edge-on galaxies, and for a trend of increasing flaring with Hubble type (in the sense that 
earlier types show more flaring). For SO galaxies (de Vaucouleurs type ~ —2), they predict 
a change in scale height with radius of ~ 9 — 10%. For Sc galaxies (de Vaucouleurs type ~ 
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5), they predict a change in scale height with radius of ~ 1 — 4%, which is consistent with 
our result of a statistically insignificant variation between the inner and outer parts of the 
disk for NGC 4244. It is hard to understand this variation with Hubble type if the flaring 
is caused by the factors discussed above, de Grijs & Peltier fit only regions between 1.5 and 
4.0 scale lengths, so the gravitational effect of the bulge can be ignored in their study. 

For earlier type galaxies, however, the determination of disk flaring is more difficult: 
as noted by |de Grijs fc Peltier (1997)| , galaxy interactions, warps, and thick disks can also 



cause an apparent change in scale height. An incorrect bulge-disk decomposition will have 
the same effect. To distinguish flaring from the presence of a thick disk requires fainter 
surface photometry than a survey project like de Grijs & Peltier's can achieve. Viewed 
from outside, the Milky Way's thick disk would be visible at /i/ = 23 magnitude arcsec"^ 
( [Morrison 1999| , assuming i? — / = 0.5), which is the surface brightness level at which 



de Grijs & Peltier's single disk vertical profiles (their Figure 3) begin to deviate from a 
single exponential fit. 



As discussed in Section pT6| , the adopted inclination, which can be an indication of a 
warp, can have a large effect on apparent fiaring. Kinematic inclination determinations, 
such as Olling's for NGC 4244, do not exist for many galaxies, making the detection of some 
warps and correct modeling of the disks more difficult, de Grijs & Peltier do not address 
the presence of warps in the galaxies they observed. 

In conclusion, the question of how the thickness of stellar disks vary with radius is a 
thorny one to approach observationally, especially in early type galaxies where the existence 
of additional components like bulges and thick disks makes accurate separation of the thin 
disk light very difficult (see Morrison et al. 1997 for a discussion of how the problem is 
exacerbated if one- dimensional profiles are used instead of two-dimensional fits). Our data 
show a statistically insignificant increase in scale height between the inner and outer parts 
of the disk of NGC 4244, and it is possible that this is partially due to the small warp. This 
highlights the need for more theoretical work on disk heating mechanisms. 



4.4. Is there a second component? 

A single exponential disk (Equation |]) is a good fit to the data down to 27.5 magnitude 
arcsec"^ (see Figure |^). We do not see evidence of an extra component. However, our data 
may not rule out a faint component resembling the Milky Way's halo, which has a very 
centrally concentrated space distribution, with density p varying as r~^ or r~^'^. NGC 4244 
has a nucleus with total luminosity 2% of its disk's, which we calculate using the disk 
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central surface brightness and bulge effective intensity given in Baggett et al. 1998, and a 
bulge effective radius of 4" (71 pc)Q. What if this nucleus was the central part of such a halo? 
We assume that the core radius of this putative halo is two pixels (71 pc)[], and calculate 
what brightness a power-law halo which fit smoothly onto the luminosity distribution of 
the nucleus would be. On the minor axis, 2 kpc above the plane, the r~^ halo would have a 
brightness /ir = 26.5 magnitude arcsec"^, which our data rule out, but the r~^^^ halo would 
have a brightness /ir = 28.4 magnitude arcsec"^, below our detection limit. 

It is interesting to compare NGC 4244 with its Local Group analog, M 33, which has a 
bulge-to-disk ratio of ~0.02 ( [Bothun 1992| ). While it is difficult to search for halo field stars 
in M 33 because of its nearly face-on inclination and badly warped disk, there is evidence 
for a small number of globular clusters with halo kinematics. [Schommer et al. (1991)| show 



that M 33 has ~20 old clusters in a kinematically hot configuration (see also |Sarajedim 



et al. 19981) . Whether these clusters are better associated with M 33's disk or tiny bulge, 
their kinematics suggest that any field stars associated with them would have a larger scale 
height than the thin disk. 

There have been no searches for globular clusters in NGC 4244, but we see no evidence 
for an extended field star population in this galaxy. A globular cluster search in NGC 4244 
would help us decide whether it truly has no halo. 

NGC 5907, a more luminous Sc, has a second component, which does not resemble 
either the Milky Way's halo or thick disk QSackett et al. 1994| , Morrison et al. 1994). 



5. CONCLUSIONS 

We have obtained deep R band surface photometry of the small, edge-on Scd galaxy 
NGC 4244. 

1. We find that NGC 4244 is well-described by a single exponential disk with a scale 
height of 246 ± 2 pc, a scale length of 1.84 ± 0.02 kpc, and a radial cut-off at 5.4 disk 
scale lengths. The thin disk parameters, h^, hz, and Rmax for NGC 4244 are smaller 
than those of NGC 5907 by about 50% (cf. Morrison et al. 1994). 



^Baggett et al. give a bulge effective radius of 400"(7.1 kpc). We assumed that this was actually 4"(71 
pc). 

^the core radius of the Milky Way's halo is not well determined but is probably less than 500 pc based 
on counts of RR Lyraes in bulge fields 
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2. NGC 4244 is a good galaxy for investigating stellar disk flaring because it does not 
have a significant thick disk and because it has no large companions. We find, in 
agreement with [van der Kruit fc Searle (1981a)| , statistically insignificant flaring in 



the stellar disk. The warp in NGC 4244 may exaggerate the marginal disk flaring that 
we see. 

3. We have placed stronger limits on the existence of a cutoff in the stellar disk. 

4. Unlike the more luminous Sc galaxies M 33 and NGC 5907, NGC 4244 does not show 
a significant second component. 



A. The Error Model 

It is necessary to have accurate error estimates for both the 1-D and 2-D fitting. 
Unlike earlier photographic work, deep CCD surface photometry allows us to quantify 
measurement errors. Measurement errors arise from CCD behavior such as readout noise 
and fiat-fielding, as well as from intrinsic galaxy variations. Each error contribution will be 
addressed below. 



A.l. Readout Noise 

The readout noise per exposure is 1.2 ADU. By combining 66 galaxy images with a 

median, we are able to reduce the effective readnoise to 

1 22 
i?e// = 1-2 ADU ^^ = 0.18 ADU. (Al) 

vNg 

Averaging the 66 frames would have reduced -Re// to 0.15 ADU; however, combining with 
a median is preferable for the removal of cosmic rays and other data defects. The factor 
of 1.22 is because of the reduced efficiency of a median (see Morrison, Boroson & Harding 
1994). 



A. 2. Photon Noise 



For C ADU in a given pixel, the photon noise is JC/g, where g is the gain. Combining 
66 galaxy images using a median reduces the photon noise to 

- 1-22 VC^ 

^Poisson 7^^ ^- \^^} 

V66 ^ 
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In the final combined galaxy image, the contribution to the error from photon noise is 
0.137 VCADU. 



A. 3. Flatfielding Errors 

In principle, the only limit to the precision of the combined flat field image is the 
photon noise in the individual fiat field images. This small-scale variation is 

VC^l.22 1.22 

where Cs is the number of counts in the final, combined master sky fiat image, g is the 
gain, Np is the number of individual sky fiats used to make the master sky flat, and Nq 
is the number of individual galaxy images used to make the final galaxy image. The sky 
counts Cs = 1998.0 ADU, the gain is 1.2 ADU, and the numbers of sky and galaxy images 
are 70 and 66, respectively. 

In practice, the small-scale flat fielding errors are not the only flat-fielding error 
we have. There are also large-scale variations which arise from the variation of the 
sky brightness across an image and from the wings of bright stars which have not been 
completely removed by combining the individual sky flats. To estimate the large scale 
variations o"z//, we divided our sky flats into two subsets of 12 each and constructed two 
mini-master sky flats. We then used IRAF's blkavg to bin the mini-masters into 100x100 
pixel bins, to remove small-scale variations, and divided one mini-master into the other. 
The standard deviation of the divided images was 0.085% Cs- The large-scale variation 
aiff is reduced by 1/v^ because the mini-masters were constructed from half as many sky 
images as the master and by another 1/2 because of the division of the two mini-masters. 
The large scale flat-fielding variation o"/// = 0.03% Cs- 



A. 4. Surface Brightness Fluctuations 

Another source of error arises from intrinsic variations in the galaxy. The light in each 
pixel comes from individual stars, and, thus, is subject to counting statistics. The variance 
of the brightness fluctuations is characterized by P?onry fc Schneider (1988"]| . Following their 



Equation 10 and [Morrison et al. (1994)| , we used t = 600s, d = 3.6 Mpc, M = (for the R 
band), and mi = 22.23 magnitude arcsec"^. The variance of the brightness is 3.61^, where 
'g is the number of counts from the galaxy only. 
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Seeing has the effect of lessening the apparent brightness fluctuations. We modeled the 
effect of seeing empirically, by creating an IRAF image with known aorig and smoothing 
the image by a Gaussian of similar FWHM as our typical seeing. The smoothed image was 
binned in bins of 3x3 to 50x50 pixels and a recalculated for each bin size. Without the 
effect of seeing, aorig would be reduced by a factor of l/y/n, where n is the bin size. We 
calculated the residuals Aa = a — (Jorig/y/n for each bin size. A function was fit to the Act 
residuals and this function was used to correct the surface brightness variation errors for 
seeing. 



A. 5. Sky Variations 

Our error in determining the sky level was 0.8 ADU, or R = 29.3 magnitude arcsec"^, 
obtained from the histogram shown in Figure ^ 
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Figure Captions 

Fig. 1. — The environment of NGC 4244: galaxies in the B 4 and B 5 galaxy groups, plotted 
in supergalactic coordinates, from Tully 1988. Galaxies brighter than Ms^tot = ~ 17.5 
(galaxies which are at least as bright as NGC 4244) are shown by solid symbols, less luminous 
galaxies by open symbols. Total range in luminosity is from MB,tot = —13.0 to —19.4. We 
have circled NGC 4244. 

Fig. 2. — Our final galaxy image, formed from 70 10-minute exposures. The orientation is 
N is up, E is to the left. The field shown here is 35'by 30'. 

Fig. 3. — Residuals from our photometric calibration are plotted as a function ofV~R color 
of the standard star. Since we observed in only one filter {R), we are unable to correct for a 
color term in our calibration, but the plot shows that the effect of neglecting the color term 
is less than 0.10 mag over the whole range of stellar colors. 

Fig. 4. — The galaxy image, rotated counter-clockwise by 42° so that the major axis is parallel 
to the X axis, with our mask superimposed. The left side of the figure is the southwestern 
side and the right is the northeastern side. 

Fig. 5. — Histogram of sky values from 310 large bins, located away from the galaxy. 

Fig. 6. — The stellar PSF, normalized to have the same flux as the central region of 
NGC 4244. At a radius of 2 kpc (56 pixels) the scattered light from the center of the 
galaxy contributes 28.5 magnitudes arcsec"^. 

Fig. 7. — An example of a contour plot of starting x^ values for fixed hz = 6.6 pixels (233 
pc). Contours are logarithmically spaced, with an interval of 0.1. The final estimate for h^ 
and Ic is shown with a cross. 

Fig. 8. — Histogram of hz values estimated by our bootstrap code. The a measured from 
these data gives our error estimate for hz- 

Fig. 9. — 1-D vertical profiles of NGC 4244, with our best fit exponential disk model shown 
as a solid line. The mean sky brightness of i? = 20.9 is marked with an arrow. Adjacent 
profiles have been shifted by 1.4 kpc for clarity. Open circles are used for data to the right 
of the minor axis, filled for data to the left. 

Fig. 10. — 1-D horizontal profile of NGC 4244, at ± 760 pc from the major axis. Points to 
the left of the minor axis on Figure 4 are shown by open circles and points to the right of the 
minor axis are shown by closed circles. The point at 9.5 kpc appears lower than it should. 
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because it straddles the edge of the galaxy and so includes a section outside the disk cutoff. 

Fig. f f . — f-D horizontal profile of NGC 4244, at ± 955 pc from the major axis. All four 
quadrants of Figure 4 are plotted using different symbols: closed circles and triangles are 
left of the major axis, open to the right, triangles below the major axis and circles above. 
The warp is clearly visible by ~5 kpc from the minor axis. The full line is the model. 

Fig. 12. — 1-D vertical profiles of the outer disk, showing the disk cutoff at Rmax = 10 kpc. 
The filled points lie within Rmax and the open points outside of Rmax- By 9.9 kpc, the warp 
has a large enough effect that points at 9.9 kpc do not lie on the model line. Open circles 
are used for data to the right of the minor axis, filled for data to the left. 
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Table 1. Errors in a bin with 500 ADU 
from galaxy & 1998 ADU from sky 



Source 


ADU 


% 


readout noise 


0.03 


0.00 


Poisson stats 


0.75 


0.03 


small-scale flat fielding 


0.10 


0.00 


large-scale fiat fielding 


0.60 


0.02 


surface brightness 


4.60 


0.19 


fluctuations 






large scale sky variations 


0.80 


0.03 


+ sky determination 







Table 2. Results of the 2-D Fitting. Parameters with errors quoted were estimated, the 

rest fixed. 



Parameter 



Value 



Correlation coefficients 
hR hz /io inc. Rr, 



hR 1.84 ± 0.02 kpc 

hz 246 ± 2 pc 

/xo (observed) 20.47 ± 0.05 R magnitudes arcsec"^ 

inc. 84.°5 



1.000 



Rn 



10.0 kpc 



0.442 


-0.678 


-0.338 


-0.465 


1.000 


-0.512 


-0.520 


-0.110 




1.000 


-0.228 


0.306 






1.000 


-0.039 
1.000 



Ho (face-on) 
liQ (edge-on) 



22.07 R magnitudes arcsec ^ 
19.90 R magnitudes arcsec"^ 
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Table 3. Comparison to Previous Studies 



Author 



hz (pc) Hr (kpc) Rmax (kpc) 



jiQ (magnitude arcsec 



van der Kruit & Searle (1981a)' 



209 



1.87 



Kodaira & Yamashita (1996) 158 

OUing's (1996) reanalysis 210 

of van der Kruit & Searle (1981) 



3.15 
2 



9.91 



10 



20.6 (-Bj, edge-on, sech^) 
^ 19.9 (-Bj, edge-on, exponen 
^ 20.1 (5, edge-on) 
^ 18.6 (-R, edge-on) 

21.6 ± 0.07 {B, face-on) 
^ 17.6 ± 0.07 (i?, edge-on) 



This study 



246 ± 2 1.84 ± 0.02 



10.0 20.47 ± 0.05 (i?, observed) 

19.90 (i?, edge-on) 
22.07 (i?, face-on) 



^Revised for our adopted distance of 3.6 Mpc 



Table 4. Comparison of Thin Disk Parameters between 
NGC 4244, NGC 5907, & the Milky Way 



Galaxy 



hz (pc) 



-l^max I '^R 



NGC 4244 
NGC 5907 

Milky Way^ 



246 
470 
300 



5.4 
3.8 

5.7 



^References for Milky Way values can be found in Morrison 
et al. (1997) 
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